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We study the Quantum Field Theory of nonrelativistic bosons coupled to a 
Chern-Simons gauge field at nonzero particle density. This field theory is relevant 
to the study of anyon superconductors in which the anyons are described as bosons 
with a statistical interaction. We show that it is possible to find a mean field solution 
to the equations of motion for this system which has some of the features of bose 
condensation. The mean field solution consists of a lattice of vortices each carrying a 
single quantum of statistical magnetic flux. We speculate on the effects of the quan- 
tum corrections to this mean field solution. We argue that the mean field solution is 
only stable under quantum corrections if the Chern-Simons coefficient N = 2tt9 / g 2 
is an integer. Consequences for anyon superconductivity are presented. A simple 
explanation for the Meissner effect in this system is discussed. 
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It is now well established that a gas of charged anyons in two spatial dimensions at 
T = is a superconductor. The most compelling evidence for this is from the calculation 
of the current-current correlator in the Chern-Simons description of the anyon gas. This 
calculation, which was first done by Kalmeyer and Laughlin Jl| in the RPA approximation 
and since worked on by many authors [@-f§], finds a massless pole in the current-current 
correlator when the anyons are electrically neutral. This implies that a neutral gas of 
anyons is a superfluid and, as a consequence, a charged gas of anyons exhibits a Meissner 
effect and is a superconductor. 

In the field theoretic work described above one treats the anyons as fermions coupled to 
a Chern-Simons (statistical) gauge field. It is this coupling which gives the fermions anyonic 
statistics. It is well known that anyons can be described either as bosons or as fermions f| 
coupled to a Chern-Simons field. The coupling, in each case can be adjusted to give the 
desired statistics. Although most of the field-theoretic work treats the anyons originally as 
fermions there are many interesting approaches to anyon superconductivity which begins 
with anyons as bosons H0]. Unfortunately it has not been possible to find an appropriate 
field-theoretic analysis to treat a gas of anyons using the scalar Chern-Simons description. 

The main problem blocking such a description is that a noninter acting system of bosons 
at nonzero density undergoes bose condensation (at T = in 2 + 1 dimensions). This makes a 
mean field (and consequently an RPA) description of the system coupled to a Chern-Simons 
field very difficult. It is the goal of this paper to study this question in detail and to derive 
the consequences of a mean field analysis for anyons which are described as scalars coupled 
to a Chern-Simons field at nonzero density. We shall see that no homogeneous mean field 
solution exists and that in mean field theory the system is unstable to the formation of a 
kind of lattice of vortices |8|J|. We shall present several arguments which suggests that this 
system of vortices should be a superconductor. This approach differs from that of Refs. 
H|7] in that we seek a mean field solution to the scalar field theory which can potentially be 
improved by RPA or other methods. 

We thus consider a gas of nonrelativistic anyons with mass m which can be described 
by a complex bosonic field ip coupled to a Chern-Simons "statistical" gauge field A^. We 
introduce a chemical potential fi and assume a Lagrangian of the form 

# f {do - igA ) ^ - J-| (V - igA) - 6 -e^ x A^A x - - (VV) 2 + (1) 

Notice that besides the interaction introduced by the Chern-Simons field we have introduced 

an additional point-like interaction | uMVM ■ The reason for this is that in the absence of 
the Chern-Simons coupling, the above system is known to bose condense. The expectation 
value of ip is < ip > 2 = fi/X. Thus, as is well known, the system is unstable (unless fi = in 
which case any density is possible). The quartic term is thus a mechanism for allowing us 
to use the chemical potential as a measure of the density. We shall see however that in the 
presense of the Chern-Simons interaction this will not be necessary as the Chern-Simons 
interaction itself stablizes the system at nonzero fi. (This is even the case if the quartic 
interaction is attractive as has been pointed out by Jackiw et. al. ||). 

We begin by trying to find a mean field solution to the equations of motion derived from 
Eq. (0) 

i (d - igAo)^P = -— (V -igA) 2 i> - ^ + A (VV) V> (2) 
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gifjty = 9B = 9V x A (3) 



^- {fD$ - (AV , ) t if>) = OtijEi (4) 

where E and I? are the statistical electric and magnetic fields derived from A and Dj = 
dj — igAj is the covariant derivative. 

First note that the only convariantly constant solution to the above equations (i.e. a 
solution for which D^ip = ) is if) = 0. This is so despite the fact that the chemical potential 
term drives if) away from if) — 0. To see this we decompose if) into its amplitude p and its 
phase Q 

i/> = pe M (5) 

Dpip = now implies that d^p = and that gA^ = d^Q. It follows that B = V x A = 0. 
But Eq. (0) now forces p 2 = 0. It thus follows that ■0 = 0. The physical reason for the lack 
of a covariantly constant nonzero solution is that there is an analog of the Meissner effect 
at work. A nonzero ip which is favored by the chemical potential term forces (due to the 
Chern-Simons interaction) a nonzero magnetic field which is expelled in the presence of a 
constant nonzero ip field. 

In addition to the lack of solutions with D^ip = there are no covariantly static solutions 
with D il) = 0. D ip = now implies that d^p = and that gA = 8qQ. Using this, and the 
fact that Ei = d^Ai — diAo we find from Eq.(f|) that 

E t = d (A - d t n) oc e y (A - d&) p 2 (6) 

Since p is time independent it follows that A — Vf2 oscillates in time from which it follows 
that B = V x A depends on time and thus from Eq.(|3D p depends on time. This is a 
contradiction unless if) = 0. 

Despite the fact that no covariantly constant solution exists, it is possible to find a 
solution for which d if> = with A ^ 0. As is shown in Ref. || this can be done by 
minimizing the energy functional 



H = J d z x 
subject to the constraint (Eq.(H)) that 



(7) 



VxA= 9 -\if)\ 2 (8) 

The term proportional to Aq in Eq.(|2|) appears when Ai in Eq.(^) is varied with respect to 
if) using the constraint (§). This occurs provided Aq is given by the solution to Eq.(|4]) with 
BqAi = 0. It is now clear that a spatially constant nonzero if) (as opposed to a covariantly 
constant if) discussed above) will yield a configuration with infinite energy since the gauge 
potential for a constant magnetic field necessarily diverges at spatial infinity. 

There are now two possibilities. Either if) = is the mean field solution of lowest energy 
or there exists a lower energy non homogeneous mean field solution. We shall see that at least 



3 



for A > the latter option holds. In fact we shall explicitly construct field configurations 
with energy H < 0. 

In general it is not possible to actually find the minimum of the energy for this highly 
nonlinear system. One approach is to tune A so that there exist solutions of the Bogomolny 
[TIJ type to this system. This requires A < i.e. a repulsive interaction. Solutions of this 



type have been studied in detail by Jackiw et. al. || and by Olesen || in the absence of the 
chemical potential term. They do find exact solutions whose structure is that of a lattice of 
vortices. In this letter we discuss the more general case when an exact Bogomolny solution 
is not available. The system will be studied using variational methods. 

Our first example of a configuration with lower energy than that of the ip = solution 
consists of a single "vortex" with a total flux $. Consider, as a variational ansatz, a config- 
uration for which the radial component of the gauge potential A r = and whose azimuthal 
component 

A t ^'^ (9) 
zng r 

so that 

B = ? P 2 = (10) 
6 F 2ngr dr V ' 

Since p 2 > it follows that B is always positive and that h(r) is a monotonically increasing 
function of r. We further demand that h(r) vanish as r — > so that there is no additional 
delta function source of magnetic flux at the origin. Since the total flux is simply g J B d 2 r = 
$ x h(oo) we choose h(oo) = 1. 
As a simple example choose 

h(r) = 1 - e~^ r)0 (l + (£rf) (11) 

so that h{r) ~ (£r) 2/3 /2 as r — > 0. For a single vortex there is no need to choose $ to be 
a multiple of 27r. We shall see later that in the case of a lattice of vortices, this will be 
essential. In anticipation of this fact we shall consider here a vortex whose flux is 2n. In 
this case we see that for large r 

gA^ — > - as r — ► oo (12) 

We can eliminate the effects of this 1/r tail by letting the phase Q of the field ip (Eq.([|)) be 
given by 

n = (f) (13) 

where is the azimuthal angle. This then implies that Vfi = <p/r so that gA — Vfi falls off 
exponentially for large r. This will turn out to be essential for the multi vortex solutions. 

Using this ansatz we can now compute the energy of this configuration. For simplicity 
we consider the case for which A = 0. We shall discuss the more general case in a future 
publication. Using Eq.(|j]) the energy can be written as 
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H = Jd 2 x! [ ^ [{d lP f + p 2 (gA t - d t nf] - pp 2 } (14) 



This integral can be done numerically for arbitrary (3 and the value of (3 which minimizes 
the energy for fixed £ can be found. This value turns out to be (3 ~ 1.7. (One should not 
take this value of (3 too seriously since it applies to the single vortex case only and not to 
the case of a lattice of vortices which is of greater interest.) At this value of f3 the expression 
for H is given by 

f 

H ~ N{ — x 1.5 -u 

where the factor 1.5 is the approximate numerically obtained value. Note that for sufficiently 
small £ the energy is lowered below H = which is the value of H for if) = 0. We thus see 
that there exist lower energy solutions than the if) = solution. The energy is lowest when 
£ — > which corresponds to a very large (or diffuse) flux tube. In the case of a lattice of 
vortices we shall see that the tubes do prefer to have a definite finite size. 

It is interesting to note that these vortices differ significantly from the usual Abrikosov 
vortices found in superconductors. First of all the order parameter p vanishes outside the 
vortex. Secondly the magnetic field tends to zero as r 2f3 ~ 2 at the centre of the vortex. Both 
of these peculiar properties are a result of the Chern-Simons condition which forces B to 
be proportional to p 2 . 

We now move on to multivortex solutions which we shall show lower the energy even 
more. In fact for a uniform distribution of vortices we shall see that the mean energy density 
is less than zero. Our technique, as for the single vortex, will be to construct a variational 
ansatz for a multivortex configuration. We imagine a distribution of n vortices at locations 
r*i • • • r n . Recall from the single vortex solution that it is imperative that ip^tp vanish at 
the location of each vortex so that the integral of (<9jf2) 2 p 2 converges there. Furthermore 
no vortex can have a long range 1/r component which is uncancelled by a <9jf2 term. If 
it did, then the integral of (gAi — dVL) 2 p 2 over all the other vortices would diverge. The 
only option is for each vortex to have a flux of 2ttx an integer so that the 1/r piece can 
be cancelled at infinity. In general vortices with flux 2n (i.e. a single quantum of flux) will 
provide configurations with the lowest energy (although we shall see cases where this may 
not be so). 

A multivortex ansatz which satisfies all our requirements is 

9Mr) = x I[h(\r-r,\) (16) 

\i=l I' 'i\J j=l 

Since h — > 1 for large r we see that this ansatz behaves like a single vortex near each of the 
vortices. To cancel the long range behavior of A we choose the phase £1 of ip to be 

fi(f)=f>(r>!) (17) 
i=i 

where <p(r, r,) is the azimuthal angle which the vector r — r*j makes with some arbitrary axis. 
As a result the important combination gA — Vf2 goes like 
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g A-vn= Et^t n^-^D-i ( 18 ) 



It is now clear that the contribution of this term due to any given vortex 'C at the location 
of any other vortex 'D' vanishes exponentially as the distance between 'C and 'D' increases. 

We could now proceed by choosing a specific form for h(r) as might be given, for example, 
by Eq.([ll]) and minimizing the energy functional from Eq. flli]) . We can learn much more, 
however, by doing a general analysis of the form of the answer. Let us, for definiteness, 
imagine that the vortices are arranged on a square lattice. Let d be the separation between 
the vortices and let l/£ be a measure of the size of the vortex (as is given, for example, in 
Eq . fllT]) ) . The energy in Eq.(|T4]) has two pieces. The first piece proportional to 1/m and 
the second proportional to /i. 

Using the fact that gp 2 = 9B and the definition of N = 2ir8/g 2 and recalling that the 
total flux per vortex is 2tt we can compute the second term in Eq.flUp as 

J d 2 x pp 2 = Nnp (19) 

where n is the total number of vortices 

n = I (20) 

where V is the volume of the system. Evaluation of the first term in Eq.(|l4]) is somewhat 
more involved. The terms proportional to 1/m are functions of only £ and d and have the 
dimensions [mass] 2 . Thus this first term can be written as 

^iV< 2 x F(Zd) (21) 
where the function F is unknown until the integrals have been done explicitly. But nonethe- 



less it is a function of the dimensionless product of £ and d ||11|| . Combining these two terms 
we can write the total energy 

The behavior of x 2 F(x) for small and large x is easy to determine. If £ — > oo the vortices 
are very small relative to their spacing and we obtain simply n copies of the single vortex 
solution. As a result F — > constant and thus x 2 F(x) — > oo as x —>■ oo. On the other hand 
when £ —>■ 0, F diverges very severely since the magnetic field is now nearly constant and 
the energy thus grows like V 2 . Thus not only F but also x 2 F(x) diverges as x — * 0. As a 
result, x 2 F(x) reaches some minimum value as x is varied. Let 7] be this minimum value of 
x 2 F(x). It is clear that rj will be a number of order 1. The precise value depends on the 
precise form of h(r) and it can be computed. It follows that for this minimum value of £d 

VN ( 77 1 
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Note that d is still a variational parameter. The minimum value of H min occurs when 
d 2 = rj/mfj, in which case 



m/i 



2 N 



(24) 



mm 



-V 



2 V 



Note that the condensate provides us with a negative energy density which is thus lower 
than the energy for the ip = solution. 

We conclude that in the context of the variational calculation it is evident that a "con- 
densate" of vortices provides a good candidate for the Mean Field ground state of this 
system. 

But this is not the whole story! Let us now compute the charge Q of one of these vortices 
in the "lattice". 



There are N anyons in this vortex with a total charge Ng. This is fine if N is an integer. 
But if N is not an integer then Q is not a multiple of g and the number of particles in 
the vortex is not an integer. Although the mean field solution still exists in this case we 
should expect to run into serious trouble when we consider quantum fluctuations about this 
vortex configuration. To see this let us transport one vortex around another vortex. The 
Bohm-Aharonov phase picked up is now 



Thus if iV is not an integer these vortices are themselves anyons. If N is rational (N — p/q 
with p, q integers) then it may be possible to consider a mean field solution with flux 2nq per 
vortex. But if N is irrational then it seems clear that the quantum corrections will lead to 
infinities. It is also likely that quantum corrections will force the number of particles in the 
vortex to be an integer. They may also ruin the lattice structure of the mean field solution 
and cause the system to form a liquid. 

Let us now concentrate on the case in which N is an integer. We must consider separately 
the case when N even and when N is odd. In the case N even (which includes the case 
of semions for which N = 2) the vortices are bosons. When one vortex is taken around 
another, the phase is exp(47ii) which leads to a phase exp(27ii) under interchange. These 
bosons have a repulsive core. It is well known that bosons with a repulsive core in two 
space dimensions form a superfluid. (At nonzero temperature there is no long range order 
in 2+1 dimensions but the superfluidity persists up to some critical (Kosterlitz-Thouless) 
transition temperature.) If N is odd the vortices are fermions. It is likely that when quantum 
effects are included the energy will be lowered by combining two vortices into a single vortex 
with flux An. This new mean field vortex will then be a boson which will allow for bose 
condensation and superfluidity. 

It is possible to see from the mean field theory directly how the Meissner effect arises in 
this theory. To this end let us couple our system to a "real" fixed background electromagnetic 
field A/x- If Ao — then the mean field problem can be formulated as a search for the 
minimum of 




(25) 



igN J Bd 2 x _ 2mN 



y£ 1 if N is not an integer 



(26) 
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H= d 2 x 



1 2 

— \(di - igAi - ieAi) if>\ - £#r 



(27) 



The mean field solution which we found for e = will certainly not admit a constant magnetic 
field b = V x A. If we choose, for example, a vector potential A2 = bx\ with A\ = then 
Eq.(p7|) will have a term / d 2 xp 2 [f)x\ which clearly diverges due to the presence of the 
infinite lattice of vortices. 

The above argument is however too crude since it is possible to make new kinds of 
vortices in the presence of an externally applied magnetic field. One can simply force the 
sum of the flux of the statistical and the real gauge fields to be 2n so that 

J d 2 x [gV x A + eV x A] = 2ir (28) 

while p 2 is proportional to only the statistical magnetic field B = V x A. This will yield 
a perfectly sensible mean field solution for arbitrarily small applied magnetic field. The 
problem with this solution is that the charge of the vortex will no longer be an integral 
multiple of g 

Q= [ d 2 xp 2 = 7— / d 2 xV xA=£gN (29) 

and, more importantly, the total number of particles in each vortex will no longer be an 
integer. As discussed previously the quantum corrections should cause such a configuration 
to have an infinite energy. As a result, a constant magnetic field is expelled from this system. 

Although the bosonic picture of anyonic superconductivity looks on the surface to be 
dramatically different from the fermionic description there are many similarities. In fact 
many of the predictions of this model are very similar to those for the fermionic descrip- 
tion. For example a crude calculation of the real magnetic field generated by these anyons 
which should be observable in muon spin rotation experiments agrees quite well with the 
result obtained in the fermion analysis. Details of this calculation will be given in a future 
publication. The one major difference between this analysis and the fermionic RPA analy- 
sis is in the behavior of the system at nonzero temperature. In the fermionic analysis the 
RPA calculation alone leads to a loss of superconductivity at nonzero temperature H. In 
the bosonic case we have a conventional Kosterlitz-Thouless kind of superconductor which 
continues to be a superconductor up to some critical temperature. 

In this paper we have discussed a nonrelativistic theory of Bosons coupled to a Chern- 
Simons gauge field. It is also possible to carry out the analysis of a relativistic Chern- 
Simons model at nonzero density. The details of that analysis differ considerably from the 
nonrelativistic analysis but the conclusions remain the same. The system forms a lattice of 
vortices which expected to be a superfluid. Details of the relativistic analysis will also be 
presented in a future publication. 
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